Two interacting quasiparticles above the Fermi sea 
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We study numerically the interaction and disorder effects for two quasiparticles in two and three 
dimensions. The dependence of the interaction-induced Breit-Wigner width F on the excitation 
energy e above the Fermi level, the disorder strength and the system size is determined. A regime 
is found where F is practically independent of e. The results allow to estimate the two quasiparticle 
mobility edge. 

PACS numbers: 72.15.Rn, 71.30+h 



Recently the combined effect of interaction and disor- 
der has been studied by different groups for two particles 
in a random potential [0-^. These researches showed 
an interaction-induced enhancement of the two-particle 
localization length Ic compared to the non-interacting 
length li. For low-dimensional systems (Id, 2d) the two 
interacting particles (TIP) can propagate coherently over 
a large distance h, but still remain localized. In 3d, 
the interaction can lead to TIP delocalization in an oth- 
erwise completely localized regime. In physical systems 
however one should study the interaction effect at finite 
particle density. This type of problem is much more dif- 
ficult for both the analytical and numerical approaches. 
Up to now the only theoretical treatment of this case has 
been done by Imry . He took into account the effect of 
the Fermi sea on the interaction-induced delocalization 
of two quasiparticles. According to Imry's estimate in 3d 
the two quasiparticle mobility edge is lower compared to 
the non-interacting one. So far only the Id case has been 
treated numerically p. However, this case is of lesser 
importance since the enhancement is recovered quite far 
above the Fermi level contrarily to the 2d and 3d cases. 
Therefore it is very important to study the problem in 
higher dimensions. 

Exact numerical computations at finite density are 
quite heavy and due to that only small system sizes are 
accessible 0,1). Even though this approach has led to 
a number of interesting results, it seems that size re- 
strictions would not allow to check the Imry estimate 
since it requires a relatively large one-particle localiza- 
tion length li. Therefore we chosed another approach 
which is based on the computation of the interaction- 
induced Breit-Wigner width F of the local density of 
states of two interacting quasiparticles (TIQ) above the 
Fermi sea. This width plays an important role since it is 
directly related to the enhancement factor for the local- 
ization length K ^ Fpc7 where pc is the density of states 
coupled by the interaction . Also F strongly affects 



localization in the Id case To facilitate the numer- 
ical simulations we used the approximation proposed in 
based on the consideration of only two quasiparti- 
cles above the Fermi energy Ep neglecting all TIQ tran- 
sitions involving hole excitations below Ep. With such 
an approximation the quasiparticle lifetime becomes in- 
finite or in other words the inelastic processes are sup- 
pressed (see discussion below). In this context we are 
able to study the TIQ problem in blocks of linear size up 
to L = 30 in 2d and L = 10 in 3d which are significantly 
larger than in exact diagonalization of multiparticle prob- 
lems. Our approach can also give a better understanding 
of the problem of quasiparticle interactions in a quan- 
tum dot which has been recently addressed experimen- 
tally [|l2| and theoretically ||l|l. 

For numerical studies we chose the TIQ model with 
on-site interaction of strength U on the 2d/3d Anderson 
lattice with intersite hopping V and diagonal disorder 
homogeneously distributed in the interval [~W, W]. The 
eigenvalue equation expressed in the basis of two-particle 
unperturbed eigenstates reads 
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The Fermi sea is intro- 



the S]2(i?)-statistics on energy scale E > T 10 1. Such 
F-based approach even though not direct is much more 
efficient and allowed to get a better understanding of TIP 



where R determines the transformation between the lat- 
tice sites basis \n) and the one-particle eigenbasis (j)m 
with eigenenergies Em so that |n) = J2m ^n,m4>m- Ac- 
cordingly Xmi,m2 a-re eigenfunctions of the TIQ prob- 
lem in one-particle eigenbasis. The matrix of transi- 
tions produced by the interaction is 

71 ^7i,mi-^n,m2 rijTTi'j^ 71,7? ^ 

duced by restricting the sum in (1) to rrii ^ with unper- 
turbed energies E / > Ep. The value of Ep is deter- 
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mined by the filling factor /i which was fixed at /i = 1/4 
in 2d and [i — 1/3 in 3d. However the results are not sen- 
sitive to this choice. At small disorder this gives approx- 
imately Ep w —1.41/ and Ep w —V respectively. Due to 
the on-site nature of the interaction, only symmetric con- 
figurations were considered. By direct diagonalization of 
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the model (1) we computed the local density of states 



Pw{E ~ Emi — Ern2 ) 
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This function characterizes the probability contribution 
of eigenfunction Xml,m2 with eigenenergy E\ in the un- 
perturbed basis \(j)mi4'rn2)- Generally we found that 
it has the well-known Breit-Wigner shape pbw{E) = 
T/{2tt{E'^ + rV4)) (see Fig.l) where the width T de- 
pends on the parameters of the model. Our main aim 
was to investigate this dependence on the system size, 
the interaction strength and the TIQ excitation energy 
e = E — 2Ep above the Fermi sea. 
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FIG. 1. Energy dependence of the rescaled Breit-Wigner 
width r/Fo in 2d. Direct diagonalization (DD) data at 
W/V = 1: U/V = 0.6 with L = 15 (triangles), L = 20 
(squares); U/V = 1.5 and L = 20 (diamonds). Fermi 
golden rule (FGR) data: W/V = 1 with L = 20 {+), 
L = 25 (x); W/V = 0.5 with L = 15 (*). The straight 
line r(e)/ro = Ce/V with C = 0.52 shows the Imry esti- 
mate. Upper insert: the same in a log-log scale with FGR 
data at higher disorder W/V = 3 (left) and W/V = 5 
(right) triangles (L — 30). Lower insert: pw vs. E for 
L = 20,W = V = 1,U = 0.6, e = 0.4 fitted by pBw with 
r = O.lSFo (solid curve). 

The results for the 2d, 3d cases in the regime of weak 
disorder are presented in Figs.l, 2. For sufficiently high 
excitation energy e, the restriction imposed by the Fermi 
sea becomes unimportant and the width F(e) tends to 
the value Fq = C/^/ {VL'^) which corresponds to the re- 
sult obtained with ergodic eigenfunctions ||l|,||j9). In this 
approach the transition matrix elements have a typical 
value J7g — U^Q^ ~ U'^/L^'^ and the transition rate is 
given by the Fermi golden rule with Fq ^ U^Q^pc and 
Pc ~ L'^'^/V. The presence of the Fermi sea modifies 
this density which becomes energy-dependent Pc(e) ~ 
L^'^e/V^. As a result the width F drops with decreas- 
ing energy as |2| 
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enough disorder as long as L < li. Hence F is indepen- 
dent on the disorder strength W. Indeed for li ^ L 
this estimate is in good agreement with the numerical 
data presented in Figs.l, 2 with C = 0.52 (2d) and 
C = 0.3 (3d). Most of the data for F in Figs.l, 2 were 
obtained by direct diagonalization of the model (1). An- 
other way to determine F without computation of the 
TIQ eigenstates is based on the Fermi golden rule which 
should remain valid for moderate interaction strength. 
This approach gives F = "^tt^,^,^^,^,^ \ UQ,ni,m2;m[.m'2 P 
6{e + 2Ep — Ejn>^ — E^i^) in terms of the transition ma- 
trix elements between one-particle eigenstates only. Here 
e = Efni +Em2 ~2Ep and to improve the statistics we av- 
eraged over different mi. 2 with approximately the same 
e. As can be seen in Figs.l, 2, both methods are in 
good agreement for interaction strength U < 1.5V. An- 
other confirmation of the validity of the golden rule is 
the [/^-dependence of F obtained by direct diagonaliza- 
tion (Fig. 2). Both approaches confirm also the scaling 
F oc L^'' valid for weak disorder. We used up to 100 re- 
alizations of disorder for the Fermi golden rule approach 
and up to 500 for direct diagonalization. 




This behaviour was assumed to remain valid for weak 



FIG. 2. Same as Fig.l in 3d. DD data at U/V = 1.2 and 
W/V = 2 with L = 4 (o), L = 5 (squares), L = 6 (dia- 
monds) and L = 7 (triangles). FGR data: W/V = 2 with 
L = 10 (x); W/V = 1 with L = 8 (+). Here C = 0.3. 
Upper insert: the same in a log- log scale with FGR data at 
higher disorder W/V = 6 (left) and W/V = 10 (right) tri- 
angles (L=10). Lower insert: FGR data for F vs. U/V at 
W/V = 2, L = 6, F = e = 0.5 (solid line: F = 0.3ro). 

The situation becomes more intricate at higher dis- 
order. Here our results show that F becomes much less 
sensitive to e variation (Figs. 3). There is a clear tendency 
that at still moderate disorder F becomes practically in- 
dependent on e which has been varied over one order of 
magnitude. In the 3d case, such behaviour takes place 
even in the delocalized regime W < Wc ~ 8.2V. The 
data even indicate a small growth of F(e) with decreas- 
ing e at VF > 6V. At high disorder li decreases and be- 
comes comparable or even less than L. In this situation 
the eigenstates are no longer ergodic in the block and the 
scaling F oc L~'^ is no more valid. In the limit 1 < li <^ L, 
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it is natual to expect another scaling F cx . To check 
this scahng we computed the inverse participation ratio 
(IPR) C ~ ^1 which allowed to calculate the ergodic value 
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FIG. 3. Energy dependence of the rescaled Breit-Wigner 
width F/Fi in 2d (left) and 3d (right). FGR data in 2d; 
W/V = 1 with L = 20 (o) and L = 30 (points); W/V = 2 
with L = 20 (open squares) and L = 30 (full squares); 
W/V = 3 with L = 20 (open diamonds) and L = 30 (fuU 
diamonds). DD data in 2d at L = 20, U /V = 0.6: W/V = 1 
(+), W/V = 2 (x) and W/V = 3 (*). FGR data in 3d at 
L = 10: W/V = 2 (o), W/V = 6 (squares), W/V = 10 (dia- 
monds) and W/V — 14 (triangles). DD data in 3d at L = 6, 
U/V = 1.2: W/V = 2 (+) and W/V = 6 (x). 

At sufficiently high excitation energy, the real width 
should be F '--^ Fi which gives the correct scaling with 
system size in the localized regime according to the nu- 
merical data. This would explain why in the block of 
fixed size F increases with increasing disorder (see inserts 
in Figs.l, 2). While such estimate gives the correct value 
of F at high energies e ~ 2V (Fig. 3), it however does 
not explain the change of energy dependence with dis- 
order. We should note that even in this unusual regime 
both the direct diagonalization and the Fermi golden rule 
computations give the same values of the width F. 

It is clear that the change of the energy dependence 
of F cannot be explained in the framework of ergodic 
transition matrix elements Ug ^ {U/V)'^A^/V where 
A ~ V/L'^ is the one-particle level spacing in the block 
of size L. At the same time the numerical results for the 
two-particle density of states p2 definitely show that it 
increases linearly with the excitation energy e as p2 ^ 
e/A^. Therefore the only possibility is that at higher 
disorder the ergodic estimate for Us is no longer valid. 
Indeed from the theory of quasiparticle lifetime in dis- 
ordered metals and quantum dots p^p"5HT^ it is known 
that the diffusive nature of the dynamics should be taken 
into account. For excitation energy much bigger 

than the Thouless energy E^., the quasiparticle decay rate 
is Yd ^ U^P3 A{Ue/VEc)'^^'^ where pa - P2e/A is the 
density of three-particle states composed of two particles 
and one hole in the final state. In the other regime rele- 



vant for the metallic quantum dot A < e <C -Be this rate 
is Yd A{Ue/VEc)'^ This shows that the matrix 

elements ^ Yjj / p^ are not always given by the ergodic 
estimate in agreement with recent results | [T^ . The differ- 
ent nonergodic regimes can be described by the following 
approximate expression nTi 
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where g = Ec/A is the conductance assumed to be 
much bigger than one. According to (4) the TIQ width 
F ~ UsP2 increases with disorder W even in the metallic 
regime since Ec = D/L'^ ^ V^ /{WL)"^ with D being the 
diffusion constant. The ergodic estimate for U"^ is recov- 
ered for g > Ep/A ^ V/A corresponding to very 
weak disorder. While the exact numerical coefficients 
in (4) are not known, it gives the energy-dependence 
is in agreement with the numerical 
data for c? = 2 (Fig. 3) but in 3d the data indicate an 
algebraic dependence with power a < (a ~ —0.2 for 
W = 6V and a ^ -0.3 for W = lOV) instead of the 
theoretical value a = 1/2. There can be different reasons 
for this discrepancy. One case W = lOV^ corresponds to 
the localized regime while the theory requires a metallic 
behaviour. The other case W — 6V, even though still de- 
localized, is quite close to the critical value Wc- Our data 
indicate that in the metallic regime with 1 < W/V < 6 
the power a smoothly changes from 1 to -0.3. 

Surprisingly at present there are no theoretical pre- 
dictions for C/^ not only near the critical value Wc but 
also in the localized regime with /i 3> 1. It seems nat- 
ural to make the assumption that in the localized case, 
the transition matrix elements will be given by an equa- 
tion similar to (4) with (7 « 1 since in a block of size 
li the Thouless energy is i?c ~ A ~ V/lf. This gives 
F ~ Fi(e/A)" where a has replaced the theoretical value 
d/2 — 1 valid in the metallic regime. We will assume that 
for (i > 2 the exponent | a |< 1. In 3d the TIQ mobility 
edge em2 is defined by the condition k = Ypc > 1 
Since pc ^ e/A^ the above expressions for F give 
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where Cmi ~ Vl^^^^ is the one-particle mobility edge. 
The one-particle critical exponent is ^ « 1.5. Due to 
that for J7 ~ y, the edge em2 ^ emi- The above result 
(5) gives for a much smaller value than the one given 
by the Imry estimate [Q. The main reason for this is 
that the transition matrix elements in the block of size l\ 
where 5 « 1 are much larger than their ergodic value used 
in 1^ . The condition that the TIQ delocalization border 
in C/ at e - F is the same as for TIP ([/ > V/l^l"^ 
gives a = 0. 

The numerical results for the dependence of k = Ypc on 
e are presented in Fig. 5. To determine numerically the 
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density of coupled states pc ^ elf'^/V'^ in the localized 
regime we computed it taking into account only those 
TIQ states which give contributions larger than 30% of 
the value of F (the data were not sensitive to the cut-off 
value). The density pc defined in this way is independent 
of the system size when L > li. The data show that for 
W = 14:V there is no TIQ delocalization {k < 1). How- 
ever closer to the one-particle delocalization border but 
still above it {W > Wc) the value of k becomes bigger 
than one at moderate excitation energies and TIQ delo- 
calization should take place. Further investigations are 
required to check more accurately the theoretical predic- 
tion (5) for the mobility edge em2- 




FIG. 4. Energy dependence of the TIQ delocalization pa- 
rameter K = Tpc in 3d (FGR data) at U/V = 1.2: L = 7 
for W/V = 14 (open trianges); L = 10 for W/V = 2 (x), 
W/V = 6 (points), W/V = 8.4 (squares), W/V = 10 (dia- 
monds) and W/V = 14 (full triagles). The dashed line shows 
the TIQ mobility edge k — Fpc = 1. 

In our model (1) all inelastic processes involving hole 
excitations have been suppressed so that F^j = 0. In 
a general case when no special physical effects suppress 
F/5 the ratio Fd/F ~ P3/P2 e/A is larger than one 
for e ^ A. This means that the quasiparticle lifetime is 
shorter than the typical time of TIQ transitions. However 
according to recent predictions ||l^ the value of F^ in a 
quantum dot is zero for excitation energies e < ei ^ 
for U V. This result was derived for metallic quantum 
dots with 5 3> 1. It is not clear what will be the re- 
sult for dots with g ^ 1 (the blocks of size li should be 
considered for TIQ delocalization |^) but it is natural to 
assume that ci ^ A ^ V/lf. This border is parametri- 
cally comparable to the TIQ mobility edge em2 (5). The 
numerical value of the ratio tm2/^i is not known. It is 
possible that it is less than one due to the 1/3-drop of 
energy per particle in the one-particle decay process. In 
this situation the quasiparticle interaction can be more 
important than the one-particle decay. 

In conclusion we analyzed the interaction-induced de- 



localization of two quasiparticles above the Fermi sea 
in the approximation where inelastic processes are sup- 
pressed. According to our numerical results and semi- 
analytical estimates, the TIQ mobility edge is below the 
one-particle edge (€m2 < ^mi)- The obtained value of 
e„i2 is even much smaller than the value proposed in Q| 
due to the fact that transition matrix elements are much 
bigger than their ergodic value near the Fermi level. 
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